Mon. Not. R. Astron. Soc. OOO.ITlfTTIO Printed 20 May 2011 (MN I*TeX style file v2.2) 



Detailed cluster lensing profiles at large radii and the 
impact on cluster weak lensing studies 



Masamune Oguri and Takashi Hamana 



O 



00 



o 

U 
o 



> 
o 

\o 
o 



Division of Theoretical Astronomy, National Astronomieal Observatory of Japan, 2-21-1 Osawa, Mitaka, Tokyo 181-8588, Japan. 



20 May 2011 



ABSTRACT 

Using a large set of ray-tracing in iV-body simulations, we examine lensing profiles 
around massive dark haloes in detail, with a particular emphasis on the profile at 
around the virial radii. We compare radial convergence profiles, which are measured 
accurately in the ray-tracing simulations by stacking many dark haloes, with our 
simple analytic model predictions. Our analytic models consist of a main halo, which 
is modelled by the Navarro-Frenk- White (NFW) density profile with three different 
forms of the truncation, plus the correlated matter (2-halo term) around the main 
halo. We find that the smoothly truncated NFW profile best reproduces the simulated 
lensing profiles, out to more than 10 times the virial radius. We then use this analytic 
model to investigate potential biases in cluster weak lensing studies in which a single, 
untruncated NFW component is usually assumed in interpreting observed signals. We 
find that cluster masses, inferred by fitting reduced tangential shear profiles with the 
NFW profile, tend to be underestimated by ~ 5 — 10% if fitting is performed out to 
^ 10' — 30'. In contrast, the concentration parameter is overestimated typically by 
^ 20% for the same fitting range. We also investigate biases in computing the signal- 
to- noise ratio of weak lensing mass peaks, finding them to be < 4% for significant mass 
peaks. In the Appendices, we provide useful formulae for the smoothly truncated NFW 
profile. 
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1 INTRODUCTION 

, Cosmological applications of clusters of galaxies rely on the 
■ inference of cluster masses. A challenge lies in the fact that 
clusters are dominated by dark matter, which accounts for 
, ~ 80% of the mass of the universe. While the distribution 
of dark matter in clusters can accurately be predicted using 
A'^-body simulations, it is difficult to relate the dark matter 
distribution with observable characteristics of clusters, such 
as gas temperatures, optical/X-ray luminosities, and mem- 
ber galaxy distributions, because it involves many compli- 
cated physical processes associated with gas cooling/heating 
and star formations. Indeed, recent cosmological constraints 
from X-ray analysis of clusters has mainl y been limited 
by various astrophysical uncertainties (e.g.. iVikhlinin et al.l 
l2009l : lMantz et al.ll201ol ). 

Weak gravitational lensing provides a direct probe 
of dark matter distributions in clusters. This method 
makes use of coherent tangential distortion of distant back- 
ground galaxies induced by deep potential wells of clus- 
ters. For massive clusters, we can detect weak lensing sig- 
nals out to virial radii, which allows direct measurements 



of virial masses (e.g., iDahle' '2006': 'Broadhurst et all l2008l : 
[Mahdavi et al.| 12003: [Okabe fc Umctsu 200^; .Qguri et all 



53; 

d), 



'" 20091: lUmetsu et al.l |2009| . l201ll : iMedezinski et al.l l2010l : 
lOkabe et al.ll2010l l. Even for less massive clusters, one can 
measure an average mass of a given cluster sample by 
adding up lensing signals for clusters in a given sam- 
ple. This stacked lensing technique has been applied to 
the Sloan Digital Sky Surve y data to obtain constraints 
on cosmol ogical parameters ("Mandelbaum & Sell a^ '200 
IJohnston et al. 2007; Sheldon ct al. 2009; Rozo et al. 201 
and will offer a promising mass-calibration method in future 
wide-field optical imaging s urveys d Oguri fc Takad3 l201ll : 
IRozo. Wu. fc Schmidd I2OI1I ) . In addition, stacked lensing 
signals much beyond virial radii probe clustering of mas- 
sive haloes, which contains co mplementary information on 
halo masses (jOguri fc Takadali201ll 'l. 

Weak lensing also provides a fascinating way to iden- 
tify massive haloes in a way unbiased to the physical state of 
the baryon, by searching for peaks in the mass map recon- 
structed from l ensing s hear measurements (jSchneideJ 19961 : 
vanWaerbeki |2000|: iHamana, Takada fc Yoshidal |2004 



Hennawi fc Spergell I2OO5I : iMaturi et all bOOsT iPace et al.l 
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This idea has been apphed to real data to 
demo nstrate that m assive clusters can indeed be identi- 
fied^ ("Mivazaki ct al.' 



■2002", '2007; 'Cavazzi 



Soucail 

ISchirmer et al...2007,;. Hamana ct al. 2009: Kubo ct al. 



1 20071; 



Such shear-selected cluster catalogs not only offer a unique 
opportunity to study the rela tion between mass and light 
in a statistical manner (e.g., iGeller et al.ll2010l ). but also 
may provide an alternative way to constrain cosmolog- 
ical parameters through numbe r counts of peaks (e.g., 
iKratochvil. Haiman. fc Mavll2010l '). 

The interpretation of weak lensing data is usually 
made by comparing observed signals with analytical model 
predictions. For the analytic calculations of lensing prop- 
erties, it is customary to adopt a de nsity profile pro- 
posed bv lNavarro. Frenk. fc White! l| 19971 , hereafter NFW). 
The NFW profile has widely been used to extract infor- 
mation on cluster masses from weak lensing data. How- 
ever, results of such analysis can be biased if the as- 
sumption on the NFW profile is not accurate. For in- 
stance, lensing signals are determined by all matter dis- 
tributions along the line-of-sigh t, w hich cause t he scat- 
ter and bias in mass estimates (IHoek stra 2003; Dodelson 
20041 : Ide Putter fc Whitell2005l;lMarian. Smith, fc Bernstein 



20ld: iMandelbaum et all I2OI0I : iBecker fc Kravtsovl 12011 
Hoekstra et al.ll2011^ ■ which can modify lensing signals par 



ticularly at large radii. Indeed, earlier work using A^-body 
simulations has found the signific ant contribution of corre- 



Havashi fc White! 2008!: Tavio et al. 


2008 


: Cacciato et al.! 


20091: !Hilbert & White!!2010l:!Masakil 
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If the true den- 



sity profile deviates from the NFW profile, it can also cause 
a systematic bias. Understanding such bias is clearly impor- 
tant for attempts to use clusters as a cosmological probe. 

In this paper, we investigate cluster lensing profiles in 
details using a large set of ray-tracing in A'^-body simula- 
tions. We pay particular attention to lensing profiles around 
virial radii where the origin of lensing signals should change 
from main haloes to the correlated matter around clusters. 
Thanks to the large number of ray-tracing realisations, we 
can study lensing profiles out to very large radii, typically 
several tens time virial radii of clusters. We consider an an- 
alytic model that better fit the simulated signals, which is 
then used to explore potential biases of several weak lensing 
studies originating from the assumption of the NFW profile. 

The structure of this paper is as follows. In Section[2]we 
present analytic models adopted in the paper. Detailed com- 
parisons with ray-tracing simulations are made in Section O 
Using the analytic model calibrated by the ray-tracing, we 
study potential biases in various lens studies in Section |4l 
We summarise the results in Section[S| In Appendices \K\ and 
[Bl we provide formulae which should be useful for various 
calculations with the smoothly truncated NFW profile. 



2 MODELLING HALO MASS PROFILES 
2.1 Main Halo 

In this paper, we study the radial profiles of the convergence 
and tangential shear around clusters. At small scales the 
profiles are dominated by the signals from the dark haloes 
associated with the clusters (the so-called 1-halo term). The 



current most popular model of the dark halo density profile 
is the profile proposed by NFW, which is defined by the 
following form (hereafter the NFW profile): 



Pnpw('") 



(r/r,)(l + r/r,)2 



(1) 



The density parameter pa is related to the virial mass Mvir 
defined such that the average density within the virial ra- 
dius becomes equal to the nonlinear overdensity Avir, which 
we compute using the spherical collapse model (see, e.g., 
!Nakamura fc Suto|[l997l ). times the mean matter density of 
the universe. Specifically, ps is described as 



Avir(2)pm(2)Cvir 
3mnfw(Cvir) 



47rr|mnfw(Cvir 



(2) 



where Cvir is the so-called concentration parameter defined 
by 



3Mv 



47rAvir(2)Pm(2:) 

and mnfw(cvir) defined by 



1/3 



T^nivj (Cvir 



{l + x) 



;dx = ln(l + Cv 



1 + Cvir 



(3) 



(4) 



The concentration parameter is known to be correlated with 
the halo mass and redshift. When necessary, we adopt the 
following relation: 



Cvir(Mvir, z) = 7.26 



Mvir 



il + z)-"-'\ (5) 



which was derived from A^-body simulations assuming 
best-fit cosmological parameters in the Wilkinson Mi- 
crow ave Anisotropy Probe (WMAP) t hird year results 
(Mac cio. Dutton. fc van den Bosch! !2008l ). with the addi- 
tion al redshift dependence based on the simulation result 
of .Duffv et al.! (|2008 ). A well-known advantage of the NFW 
profile in lensing studies is that there are analytic expres- 
sions for the radial p rofiles of the defiection angl e, con- 
verge nce, and shear (Bartelrnann 1996; Wright & B rainerj 
!2000l ). 

However, the NFW profile is not well-defined in the 
sense that the enclosed mass diverges logarithmically. Thus 
!Takada fc jairi! (!2003a! lbl) considered lensing by the NFW 
profile truncated at the virial radius (hereafter the TJ pro- 
file): 



PT.T(r) 



(r/rs)(l + r/vs) 



re(rvir 



(6) 



with Q{x) being the Heaviside step function. The lensing 
pro perties of this profile c an also be computed analytically 
see !Takada fc Jainl!2003a! lbl). 

One potential problem of the TJ profile is that the shear 
and convergence profiles are not differentiable at the trun- 
cation radius, which causes the diverg ence in the fiexion 
profile. !Baltz. Marshall, fc Oguril |2003) proposed a differ- 
ent form of the truncation (hereafter the BMO profile): 



pBMo{r) 



Ps 



(r/r-,)(l-|-r-/r,)2 \r^ + 



(7) 



where rt is the truncation radius. In the paper we mainly 
use the following dimensionless truncation radius 
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rt 



(8) 



While both n = 1 and 2 have been considered in iBaltz et alJ 
in this paper we consider only n = 2 for simplic- 
ity. In the case, at large radii the density profile behaves as 
PBMo(f') oc r~^, and therefore the enclosed mass converges 
quickly. The three-dimensional enclosed mass of the profile 



(9) 



X (r^ + l)x {x{x + i) - t\x - i)(2 + 3a;) - 2t'^} 
+t{x + i)(r^ -I- x^) {2(3r^ - i)arctan(a;/r) 
+r{^ - 3) ln(r'(f + xf/{r^ + x^))} 



(10) 



where r = rt/vs = t^c^iy and we have fixed n = 2 in equa- 
tion. Again, the advantage of the choice of this specific 
form of the truncation is that lensing pro perties can be com- 
puted analytically (see lBaltz et al.ll2009l ). 

In equation ([7} , we adopted the same density parameter 
ps as used in the NFW profile, with Mvir as an input param- 
eter of this model. This virial mass Mvir, strictly speaking, 
is different from a true virial mass defined from the require- 
ment of the nonlinear overdensity. However, one advantage 
of this choice of the normalisation is that the "virial mass" 
of this model, Mvir, becomes rather close to the total mass 
Mtot. From the definition above, we can show that these two 
masses are related by (again for n = 2) 



Mt< 



mtot 



mnfw(Cvir) 



Mvir 



(11) 



mtot = mbn.o(oo) = ^^^a ^ i)3 

X [(3r^ - I)(7rr - - 1) + 2r^(r^ - 3) Inr] . (12) 

Typically we have Mtot /Mvir ~ 1-1.3 and M(< 
rvir)/Mvir ~ 0.9 in the parameter range we are interested 
in. 



2.2 Correlated matter 

At the large scale, the correlated matter around the cluster 
(the so-called 2- halo term) contributes to the lensing profile 
(e.g.. I Johnston et al.ll200"^ ). Using the Limber's approxima- 
tion, we compute the 2-halo term contribu tion for a cluster 
with mass M and at redshift z as (e.g., lOguri fc Takadal 
H) 



K21r(6l) 



2^''°^'''^l + z)3E„itOi(z) 



Pm{ke-z), (13) 



where kt = £/{{! + z)Da{z)}, -Da(z) is the angular diam- 
eter distance, Scrit is the (physical) critical mass density 
for lensing, Pmik) is the linear matter power spectrum, and 
hh[M) is the halo bias. We can compute the corresponding 
radial profile of the tangential shear 7t {0) simply by replac- 
ing the zero-th order Bessel function Jo{x) in equation (|13|) 



to the second order Bessel function J2{x). The power spec- 
trum Pmik) is computed usi ng the approximation o f the 
transfer function presented bv lEisenstein fc Hul (|l998t ). 

We ado pt a model of the halo bias bh{M) by 

[Tinker et al] (|2010l ). which has been calibrated using a large 
set of A'^-body simulations. The model has been presented 
as a function of several different mean overdensities A. We 
adopt a model with A = 800, a larger overdensity than the 
virial overdensity of A ~ 300, so that the details of the trun- 
cation near the virial radius do not affect the conversion of 
the masses very much. For a given set of Mvir and Cvir, we 
compute the mass for A = 800 assuming the NFW profile 
(see, e.g. jHu fc Kravtsov|[20Ci3 ). 



3 COMPARISON WITH RAY-TRACING 
SIMULATIONS 

3.1 Simulated data 

In order to study the detailed radial profiles of lensing signals 
around massive haloes, we employ a large set of ray-tracing 
simulations presented bv lSato et al.l (|2009l ). The ray-tracing 
simulations are based on 200 x 2 realisations of A'^-body 
simulations with the box sizes of 240 (particle mass nip = 
5.4x lO^°/i"^M0) and 480/i"^Mpc (m^ = 4.4x 10"/i"^Mq). 
The number of particles in each A-body run is 256^. Outputs 
of the small and large box simulations at different redshifts 
are then placed to construct a light cone of the 5° x 5° region. 
We use the standard multiple lens plane algorithm to simu- 
late lensing by intervening matter. We trace 2048^ light rays 
backward from the observer, resulting in the angular grid 
size of 0.15'. Thus our ray-tracing simulations take account 
of all density fluctuations along the line-of-sight, not just 
the matter distributions around massive haloes. While lens 
planes have been constructed out to z — 3, in this paper 
we consider only one source redshift of Zs = 0.997, which 
is a typical mean source redshift of weak lensing analysis. 
Thanks to the large number of A-body runs, 1000 inde- 
pendent ray-tracing realisations can finally be generated by 
randomly shifting the simulation boxes. The total effective 
area of 25000deg^ thus allows us to compute the average ra- 
di al profiles quite ac curately. Interested readers are referred 
to lSato et aiT j2009"l for more details of A-body simulations 
and ray-tracing technique. The assumed cosmological pa- 
rameters are the matter density ~ 0.238, baryon density 
ilb = 0.042, cosmological constant Qa = 0.762, spectral in- 
dex Us = 0.958, the normalisation of the matter fiuctuation 
(Tg = 0.76, and the Hubble parameter h — 0.732, which we 
adopt in our analytic calculations too. 

In the A'-body simulations we identify dark haloes us- 
ing the friend-of-friend (FOF) algorithm with the linking 
length of b = 0.2 to derive the total mass Mfof of each 
halo. We compute average convergence and tangential shear 
profiles in each mass and redshift bin. We consider 5 mass 
bins centred at A/fof = 10"■■^^ 10^*, 10^'*■^^ 10^''■^ and 
lO" '^'^/i"^M0, with the bin size of A(logMFOF) = 0.1, and 
7 redshift bins centred at z = 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, and 
0.7, with the bin size of Az — 0.04. The number of haloes 
used in each bin ranges from ~ 20 to ~ 20000. When stack- 
ing the profiles, we carefully choose a centre of each halo by 
identifying a peak in the convergence map in order to avoid 
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Figure 1. Bets-fit values of the dimensionless truncation radius 
Tv (eq. H]) in the BMO profile, for difi'erent redshift {upper) and 
mass (lower) bins. Filled squares plot the case when the concen- 
tration parameter is computed from equation (|5)l, whereas open 
circles are the case when the concentration parameter is fitted 
simultaneously. Dotted lines connect median values of t„ in dif- 
ferent redshift and mass bins. The median values of all bins are 
T„ = 3.0 for the fixed concentration parameter, and t„ = 2.6 for 
the fitted concentration parameter. 



any off-cent ring effect which essen tially smears out central 
profiles (see lOguri fc Takada|[2oTll ) . The average profiles are 
computed as a function of 0/0vir, where Svir is the virial ra- 
dius computed assuming Mvir = Mfof. We derive profiles 
in the radius range of 10~^ < Q jQvir < 50, with the bin size 
of log(e/evir) = 0.1. 

3.2 Fitting method 

We fit the average convergence profile in each bin using the 
models described in Sec. (2] We basically fit the average pro- 
files in the simulations with the sum of the 1-halo and 2-halo 
components, K(ff) = Kih(^) + «;2h(^), with three different 
models of the 1-halo components as described in Sec. 12.21 



The standard x method is employed for fitting, with the 
error of convergence in each radial bin derived from the scat- 
ter in the simulations. While the average profiles are derived 
for given ranges of Mfof, it is not clear whether Mfof co- 
incides with the virial mass Mvir used in the analytic mod- 
els. In fact, previous work using TV-body simulations has 



found that there are systematic b ias and cons i derable scat- 
ter b etween Mfof and Mvir (e.g., IWhitj|200ll : [Tinker et all 
l2008l ). Therefore, we perform fitting with leaving the virial 
mass Mvir as a free parameter to take account of such sys- 
tematic effect. Unless otherwise stated, we regard the con- 
centration parameter Cvir as a free parameter, although our 
main conclusions are unchanged even if we fix the concen- 
tration parameter to the value obtained from equation ([5]). 
Thus there are 2 fitting parameters (Mvir, Cvir) for the NFW 
and TJ profiles, and 3 fitting parameters (i\fvir, Cvir, t„) for 
the BMO profile. 

It has been known that very central density profiles of 
dark haloes in A'^-body simulations are not reliable because 
of several numerical effects such as the two-body relaxation 
and t he finite time step size. Ba sed on the detailed analy- 
sis of iFukushige fc Makinol (|200ll ). we estimate that density 
profiles are reliable down to ~ O.OSrvir for our A''-body simu- 
lations in which A'^ ~ 10"* particles are included in each halo 
analysed in the present paper. Indeed, our analysis results 
also indicate that convergence and tangential shear profiles 
near the very centre [Q jQvir ~ 0.01) tend to be significantly 
smaller than our analytic model calculations (see below). 
Thus, in all mass and redshift bins fitting is performed in 
the radial bins 9/9^ir ^ 0.079 to make sure that our fitting 
results are not affected by numerical artifacts. 

3.3 Result 

First we check the best-fit values of the truncation radius 
Tv (eq. [S]) in the BMO profile. Figure [T] show the best-fit 
values for different redshift and mass bins. We find the best- 
fit values to be ~ 2 — 3, showing no strong dependence on the 
mass and redshift. We consider two cases, the case that the 
concentration parameter is fixed to the value computed in 
equation ([5)l and the case that the concentration parameter 
is also fitted to the data. We find that the results of both 
cases agree reasonably well with each other, which implies 
that the best-fit values of Tv is not sensitive to how the 
concentration parameter is treated. The median values are 
r„ = 3.0 when the concentration parameter is fixed, and 
r„ = 2.6 when the concentration parameter is fitted to the 
data. 

To illustrate how well the different main halo models 
can reproduce the average profiles in ray-tracing simula- 
tions, in Figure [5] we show the comparisons of convergence 
and tangential shear profiles for a representative case. As 
clearly shown in the Figure, the BMO profile reproduce pro- 
files in the ray-tracing simulations quite successfully for a 
wide range in radii. In contrast, the NFW and TJ profiles 
are less successful in fitting the profile around the virial ra- 
dius, where the transition between 1-halo and 2-halo terms 
occurs. The NFW profile tends to overpredict the conver- 
gence profile, whereas the TJ profile clearly underpredicts 
the profile. For comparison, we fit the convergence profile 
using the 1-halo term only for the NFW and TJ profile. We 
find that the NFW profile can fit the total average profile 
quite well out to Q jQviT ~ 2, beyond which the NFW profile 
clearly underpredicts the profile. The situation is worse for 
the TJ profile for which we can see the significant discrep- 
ancy already at S/^vir ~ 1. 

The Figure also indicates that the analytic and simu- 
lated convergence profiles agree well with each other down 
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Figure 2. Example of fitting results for the mass bin Mfof = W^'^h~^ Mq and the redshift bin z = 0.4. Left panels show convergence 
profiles, whereas right panels display tangential shear profiles. Curves are best-fit results for three different main halo models, the NFW 
profile (dashed), the TJ profile {dash- dotted), and the BMO profile (solid). The average profiles from ray-tracing simulations are indicated 
by filled squares. Top: Results when only 1-halo term (i.e., no 2-halo term) is considered. The result for the BMO is similar to the NFW 
result, and therefore not shown. Middle: Results for the BMO profile. Contributions from 1-halo and 2-halo terms are shown by dotted 
lines. Bottom: Results for all the 3 profiles are compared. 



to very small radii of 9/9vir ~ 0.05, whereas the tangential 
shear profiles show large deviations already at QjQvw ~ 0.1. 
One reason for this is that shear signals are non-local. Since 
the tangential shear signal at certain radius reflects all the 
mass distributions at smaller radii, numerical effects appear 
at larger radii in the tangential shear profiles than in the 
convergence profiles. 

Figure show residuals of the convergence profile, i.e., 
the fractional difference between profiles in simulations and 
best-fit analytic models, averaged over all redshift and mass 
bins. We find that the BMO profile can fit the convergence 
profiles in simulations quite well (< 5%) for a wide range of 
radii from the core of main haloes to more than 10 times the 



virial radii. We find the NFW profile generally overpredicts 
the convergence profile by ~ 20 — 30% at S/Svir ~ 2. On the 
other hand, the TJ profile grossly underpredicts the conver- 
gence profile at Q jQviT ~ 1- Thus we confirm the result shown 
in Figure [2] that the BMO profile can best describe conver- 
gence profiles in ray-tracing simulations. Again, we check 
residuals for the case that we include only 1-halo term in 
the analytic model, and find that the NFW profile can fit 
the convergence profiles quite well out to 6/6vir ~ 1, but be- 
gins to underpredicts the profile quickly beyond the radius. 
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Figure 3. Residuals of convergence profile fitting (see Figure [2] 
for a representative example) averaged over all redshift and mass 
bins, plotted as a function of the normalised radius S/Svir- Specif- 
ically we define the residual as (Ksim — where Kgim is the 
average convergence profile from ray-tracing simulations and k is 
the convergence profile of the best-fit analytic model. Lines are 
same as Figure [2] Upper: Residuals for fitting when only 1-halo 
term is considered (see also top panels of Figure[2]l. Again, the re- 
sult for the BMO is similar to the NFW result. Lower: Residuals 
when both 1-halo and 2-halo terms are included in the analytic 
model (see also bottom panels of Figure [2]l. The thin solid line 
indicates the result when the truncation radius in the BMO pro- 
file is fixed to Tv = 2.6, the median value among fitting results 
for all redshift and mass bins (see Figure [TJ- 



Figure 4. Biases in weak lensing mass estimates from shear pro- 
file fitting, as a function of the outermost fitting radius Bout, 
derived by adopting the BMO profile as the true lensing pro- 
file. The bias is defined as [Afvir(fit) — Afvir] /A^vin where Mvir 
is the input virial mass and A/vir(fit) is the best-fit virial mass 
to the input reduced shear profile, using the NFW profile with 
no 2-halo term as a model for fitting. Note that the innermost 
fitting radius is fixed to Sin = We consider three input halo 
masses, Afvir = IO^'^H-'^Mq (solid), 3 X (dashed), 
and 10^^h~^MQ (dash-dotted). The halo redshifts are z = 0.2 
(upper) and 0.6 (lower). Thick and thin lines indicate results for 
the input truncation radii of t„ = 2.6 and t„ = 2, respectively. 
Vertical lines show virial radii for these cluster masses. 



4 IMPACT OF INACCURATE PROFILES ON 
CLUSTER WEAK LENSING STUDIES 

In most weak lensing studies, tiie NFW profile lias been 
adopted in comparing with observed lensing signals, witii- 
out including any contribution from tiie 2-iialo term. How- 
ever, tiie difference between assumed and true profiles can 
induce systematic biases in the interpretation of results. In 
this section, we investigate such systematic biases, assum- 
ing the analytic model calibrated by ray-tracing simulations 
(see Sec. [31) as a true cluster lensing profile. 

4.1 Shear profile fitting 

The most popular method to measure cluster masses from 
weak lensing data is to fit the tangential shear profile with 
analytic model predictions, for which the untruncated NFW 



profile is usually adopted as the analytic model. Here we 
estimate how cluster masses derived by such shear profile 
fitting can be biased due to the deviation of the true cluster 
lensing profile from the one computed from the NFW profile. 

In actual weak lensing analysis of observed data, we 
measure reduced shear profiles around clusters, which are 
defined by g{6) = 'yT (6) / [1 — k{6)]. For a given reduced 
shear profile, we usually derive the best-fit mass by mini- 
mizing the following x^- 

^2 ^ ^ [gjdj) - gNFw(Oi; M^ir, C^ir)f ^^^^ 

i ^ 

where g{di) is observed reduced shear profile at the radius 
9 = 9i and 5NFw(Si; Mvir, Cvir) is the corresponding analytic 
model prediction assuming the NFW profile with no 2-halo 
term. We compute g{6i) assuming the BMO profile which 
has been shown to best reproduce lensing profiles in ray- 
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Figure 5. Same as Figure (4] but biases in estimates of the con- 
centration parameter Cvir are siiown. 



ing masses are biased low by ~ 5 — 10%. They adopted 
^out = 15', 20', 25', and found a larger amount of the 
bias with increasing 6out- Thus our results are consistent 
with their results, both qualitatively and quantitatively, and 
moreover provide a physical explanation for their finding. 

In Figure [S] we show biases in concentration parame- 
ter estimates. We find that the concentration parameter is 
generally overestimated, typically by ~ 20%. Clearly this is 
related to the mass estimation bias studied above. In lensing 
analysis with the NFW profile, there is a well-known degen- 
eracy between mass and concentration parameter, and hence 
the underestimate of the mass has to be compensated by the 
increase of the concentration parameter in order to recover 
lensing signals near the cluster centre. 

A caveat here is that the input M^ir in the BMO pro- 
file differs from the true virial mass of the profile, which 
may make the interpretation of this result somewhat dif- 
ficult. In fact we can avoid this problem by considering 
larger nonlinear overdens ity A, such as A « 1600 used by 
iBecker fc Kravtsovl l|201ll ). because the enclosed masses for 
the NFW and BMO profiles adopting the same ps should 
become more similar for larger values of A. We confirm that 
our results shown in Figures |4] and [5] similarly hold for dif- 
ferent nonlinear overdensities. 



4.2 Weak lensing mass peak 

Next we discuss how the uncertainty of lensing profiles can 
affect the interpretation and statistics of weak lensing mass 
peaks. The mass peak is usually characterised by the signal- 
to-noise ratio: 

^pcak 



(15) 



where Kpoak is the filtered peak amplitude computed by 



tracing simulations (Sec. O). Assuming the uniform source 
galaxy density, we adopt the error ai oc l/6i. The x is 
computed in the range 9in < 9 < 6'out with an interval of 
A(log&) — 0.05. We fix the innermost radius to 6in = 1', 
and see how the best-fit mass Mvir differs from the input 
mass as a function of dout- In fitting, we vary both Afvir and 
Cvir to search for the best-fit parameter values. 

Figure |4] shows the results, the biases as a function of 
^out- We find that adopting a single NFW profile tends to 
underestimate virial masses for 6'out ~ 10' — 30', a typical 
outermost radii adopted in the real analysis of cluster weak 
lensing data. The exact amount of the biases appears to be 
sensitive to r„. For our canonical value of r„ = 2.6, which 
was the median value among fitting results for many redshift 
and mass bins in ray-tracing simulations (see Figure[T]) virial 
masses are underestimated by ~ 5%. On the other hand, if 
we adopt a slightly smaller value of r„ = 2, virial masses 
can be underestimated by more than 10%. If we choose very 
large values of Sout, the bias becomes smaller because of 
counteracting effects of the overpredciton of signals near the 
virial radius and the underprediction of the signals at very 
large radii where the 2-halo term contributions dominates. 
Our result suggests that we can reduce such mass estimation 
bias by restric ting the fitting rang e smal l, 9out < ^vir. 

Recently, iBecker fc Kravtsovl l|201ll l studied biases in 
weak lensing cluster mass estimates in details using a cos- 
mological A''-body simulations, and found that weak lens- 



"jpcak 



'yTi9)Q{e)2Tv9d9, 



(16) 



with Q{9) being the filter function. The noise (Jnoise is re- 
lated to the intrinsic ellipticity and the number density 
of source galaxies, for which we assume Ue = 0.35 and 
rigai ~ 30 arcmin"^. In this paper, we consider two fil- 
ter fu nctions. One is the Gaussian filter ijKaiser fc Squired 



(17) 



and t he other is a filter proposed by iHennawi fc Soergell 
(|2005l . hereafter HS05), 



1 



21og(l + a;) 



(18) 



where x = 9/9in and Xout = ^out/^in. We adopt 9in = 1' for 
the Gaussian filter, and Oin = 0.5' and ^out ~ 10' for the 
HS05 filter, which are typical values adopted in the litera- 
ture. 

We study the bias in f assuming the true profile to be 
the BMO profile. Figure |6] shows how values of v are biased 
by adopting the NFW profile with no 2-halo term. We find 
that the bias is < 4% for peaks with ~ 3 and smaller for 
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Figure 6. Contours of the biases in the signal-to-noise ratio v 
defined in equation JTS) in the mass-redshift plane. The bias is 
defined by (t'NFW ~^)/^^ where ^nfw is the signal-to-noise ratio 
computed using the NFW profile with no 2-halo term, and v is 
the true signal-to-noise ratio computed from the BMO profile 
with r„ = 2. Note that the bias is even smaller for the canonical 
truncation radius of t„ = 2.6. Dotted lines show contours for 
1/ = 3, 5, 7, and 9, from bottom to top. The top panel is the 
result for the Gaussian filter (eq. |17)). whereas the bottom panel 
is for the HS05 filter (eq. [18]). 



higher v. The bias is not so large because v is mainly de- 
termined by the lensing profiles well inside the virial radius 
where the difference of the profiles are rather small. Note 
that the result shown in Figure [6] is for the truncation ra- 
dius of Ti, = 2, which is smaller than our canonical value of 
Tv = 2.6. The bias is even smaller if we adopt Tv — 2.6. 

In some previous work, the TJ profile h as been used to 
comp ute u instead of the NFW profile (e.g.. iHamana et al.l 
|2004 ). In Figure [3 we show the bias for the TJ profile. We 
confirm that the bias has a different sign from the case for 
the NFW profile, with a similar size of < 4% for most of the 
parameter range we are interested in. 

Therefore, in either case the effect of the outer profile 
on the calculation of the peak height f is not very large, al- 
though it can be important for high-precision analytic pre- 
dictions. We note that the result presented here is an aver- 
aged effect. In practice, there is a considerable scatter on the 
amount of correlated matter around clusters, which can af- 
fect the statistics of weak lensing mass peaks in various ways. 
The detailed study of statistical properties of weak lensing 
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Figure 7. Same as Figure[6] but the biases for the TJ profile are 
shown. 

mass peaks is currently underway, which will be reported in 
a forthcoming paper (Hamana & Oguri, in preparation). 



5 SUMMARY AND DISCUSSIONS 

We have studied detailed lensing profiles and their impact 
on various cluster lensing studies. First, we compare our 
simple analytic model predictions, which consist of three 
different 1-halo terms and a 2-halo term, with a large set 
of ray-tracing in A'^-body simulations, with a particular em- 
phasis on lensing profiles at around virial radii of massive 
haloes. Our ray-tracing solves the full light propagation from 
the source plane to the observer, and thus take accounts of 
all the matter fiuctuations along the line-of-sight. We find 
that the BMO profile, which is essentially a smoothly trun- 
cated NFW profile, with the truncation radius of Ti, ~ 2 — 3 
with additional contribution from the 2-halo term best re- 
produces the lensing profiles in ray-tracing simulations. In 
contrast, the untruncated NFW profile tends to overpredict 
the convergence profiles by 20 — 30% at twice the virial radii 
when the 2-halo term is added. The TJ (hard-truncation) 
profile largely underpredicts the lensing signals at around 
the virial radii. 

Next, assuming the BMO profile with the 2-halo term as 
the "true" cluster lensing profile, we have investigated biases 
in cluster weak lens studies coming from the adoption of the 
NFW profile as an analytic model. We find that, if observed 
reduced shear profiles are fitted by a single NFW compo- 
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nent, as has commonly been done in cluster weak lensing 
analysis, we can underestimate cluster masses up to 5 — 10%. 
The bias is most significant when the outer boundary for fit- 
ting is set to ~ 10' — 30' depending on the mass and redshift 
of the cluster, or equivalently several times the virial radius. 
We can reduce the bias by restricting the fitting range small 
with the outer boundary equal or smaller than the virial ra- 
dius. On the other hand, the concentration parameter tends 
to be biased high (up to ~ 20% or even higher) to com- 
pensate the underestimate of the mass. These results based 
on the analytic mod el appear to be con sistent with recent 
numerical results by . Beck er fc Kravtsovl |2011). 

In the reduced shear fitting above, we have considered 
only intrinsic shapes of galaxies as a source of the error, i.e., 
ignored the contribution of the large-scale structure (cos- 
mic shear) to the err or fe.g.. lHoekst"rall2003l : lDodelsonll2004l : 
iHoekstra et al.l [20111 ). The error from the large-scale struc- 
ture becomes increasingly important at larger radii, and thus 
decreases the relative weight of bins at large radii, implying 
that it reduces the bias in the mass and concentration pa- 
rameter estimates. The exact amount of the reduction of the 
bias, however, depends on the source galaxy number density 
assumed. In real cluster weak lensing analysis, colour cuts 
have often been applied in order to minimise the dilution 
effect by cluster member galaxies, which significantly de- 
creases source galaxy number densities used for the analysis. 
Because of this, the contribution of the large-scale structure 
to the total error b udget tends to be subdominant in real 
data analysis (e.g., lOguri et al.ll2O10l ). but the large-scale 
structure should still contribute significantly at very large 
radu, e.g., > 20'. 

We have also examined the effect of the outer lensing 
profile on the signal-to-noise ratio i/ of weak lensing mass 
peaks. We have considered two different filters, and found 
that V can be biased by < 4% for peaks with h' >3hy adopt- 
ing a single NFW component in calculating u. Although this 
level of the bias is not so large, it can be important in de- 
tailed comparisons of observed mass peaks with theory in 
the future. 

While in this paper we have focused on biases due to 
modeling uncertainty of the outer lensing profiles, there are 
other effects that could also bias the weak lensing analysis. 
Among others, the halo triaxiality could have the most sig- 
nificant impact. Indeed, the large triaxiality of cluster-scale 
haloes suggests that the lensing signal is a strong function of 
the viewing angle (e.g.. lOguri et aLlbOOSl l. We expect that 
the effect is almost averaged ou t by proper statistical anal- 
ysis (e.g., ICorless fc Kind l2009f ). but the exact amount of 
the residual bias due to the triaxiality should be sensitive to 
how the cluster sample is selected. 

We note that the analytic model considered in the pa- 
per is subject to a potential improvement. For instance, the 
2-halo term considered in the paper is accurate in the large- 
scale limit. Near the cluster we may have to consider var- 
ious effects such as non-linearity and stochasticity of the 
halo bias for more accurate theoretical prediction. In addi- 
tion, for more robust separation of signals between the main 
halo and correlated matter, it is important to compare the 
three-dimensional matter distribution around clusters with 
projected lensing signals, which is left for future work. 

Nevertheless, our analytic models presented in the pa- 
per, which has been tested against a large set of ray-tracing 



simulations, should be useful for various aspects of cluster 
weak lensing studies as a handy model of realistic lensing 
profiles in A'^-body simulations. In particular, we expect that 
they are invaluable in the stacked weak lensing study for 
which lensin g signals can be detected much beyond t he virial 
radius ('e.g.. lMenard et al.ll201ol : lUmetsu et allbOllT ). Given 
an important c osmological informatio n coming from the 2- 
halo term (e.g.. lOguri fc Takadall201ll ). it is essential to have 
a proper analytic model which accurately predicted signals 
from the 1-halo to 2-halo terms, in order for an unbiased 
interpretation of stacked lensing data. The model is also 
useful for studying the effect of the outer profile uncertainty 
(including the contribution from correlated matter around 
clusters) on various cluster weak lensing analysis, as explic- 
itly demonstrated in the paper for shear profile fitting and 
weak lensing mass peaks. 
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APPENDIX A: EQUATIONS FOR LENSING BY 
THE BMO PROFILE 

The lens ing properties of t he BMO profile (eq. [7]) were de- 
tailed in lBaltz et al.l l|2009h . in which analytic expressions of 
lens potential, deflection angle, convergence, etc. have been 
derived. Here we reproduce key equations for both n = 1 
and n = 2 cases with the notation used in the paper. 

For n = 1, the convergence K{r) and average conver- 
gence k(< r) are given by 



K{r) 



4psr-s 



Ecrit 2(r2 + 1)2 



.{l-F(x)} 



+2F{x) 



■■L{x) 



(Al) 



k(< r) 



Scrit + l)%-2 



+ {r^ - 1) Inr + y^r'^ + \ -n + -L{x) 



,(A2) 



where x = r/rs, r = rt/rs — r^Cvir, and F{x) and L{x) are 
F{x)-- 



■ arc.tanh-v^/l — (a; < 1)^ 

(A3) 

■ arctan-v/T^ — 1 (x > 1), 



L{x) = In 



\/t'^ + x'^ + T 



(A4) 



It is also useful to show the three-dimensional enclosed mass 
'Tibmo(a;) (see eq. [S]) for n = 1. It is given by 

^2 

""^■""^"'^ " 2(r2 + 1)2(1 + s) 

X -2(t^ + l)a; + 4r(x + l)arctan(a::/r) 

-h(r" - 1)(1 + x) ln(r'(l + xf/{T^ + x"))] . (A5) 
Then the total mass rritot (see eq. [11]) is 



mtot = mbmoiooj 

„2 



-— [vrr - - H- (r^ - 1) In r] 



(r2 + 1)2 

For n = 2, K,(r) and k(< r) are given by 



(A6) 



K(r) 



Ecrit 4(r2 + 1)3 
+9,F{x) + 



2(r2 + 1) 



{l-F{x)} 



r'^ir'^ + 2:2) 



a;2 - 1 

7r[4(r2 + a;^) + + 1] 
(r2 + a;2)3/2 
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^ r^(r^-l) + (r^ + x^)(3r^-6r^-l) ^^^^ 



3(r2+x2)3/2 



, (A7) 



«(< r) = 



2{t +'^x^ -Z)F{x) 



Ecrit 2(r2 + 1)3x2 

+ i {7r(3r^ - 1) + 2r(r^ - 3) Inr} 

H ^ |-T37r(4r2 + Sr^ - 1) + 

T'^V r2 + 2:2 1^ 

+ [2r''(r2 - 3) + x\?,t* - Gr^ - 



(A8) 



The three-dimensional enclosed mass mbmo(3;) and the total 
mass mtot for n = 2 were already given in equations (|10p 
and p2p . respectively. 

From the above expressions, the tangential shear 7t(»') 
and deflection angle <?!)r(7') can be computed as 



7T(r-) = k(< r) - K(r), 
cf>r{r) = r k(< r). 



(A9) 
(AlO) 



APPENDIX B: FOURIER TRANSFORM OF 
THE BMO PROFILE 

The Fourier transform of a dark halo density profile consti- 
tute an imp ortant ingredient of the so-called halo model ap- 
proac h fe.g.. lSeliakll2000l : lMa fc Frvll2000l : [Scoccimarro et all 
I2OOI I'). In this appendix we study the Fourier transform of 
the BMO profile. 

We consider the three-dimensional Fourier transform 
of UM{r) = p(r)/Mvir, UM{k) with k being the comov- 
ing wavenumber. We note that the Fourier transform of 
lensing signals (convergence and tangential shear profiles) 
can also described by using U M(k) as a consequen ce of 
the projection-slice theorem (see lOguri fc Takadall201ll '). We 
find that um (k) for the BMO profile can analytically be ex- 
pressed as 

"""^^^ " 2m„fw(cvir)(l + r2)2a; 

2{t'' - 1)P{tx) - 2r7re-"" 



+T {vr — 231(2:)} 1 2 cos X + (r^ + l)x sin a;} 
-|-2rCi(a;) |2 sin x — (r^ + 1)2; cos x} 



(Bl) 



for n — 1, and 

UM{k) 



4mntw(Cvir)(l + r2)3a; 

2(3r* - 6t^ - l)P(ra;) - 2t{t^ - 1)xQ{tx) 

-2T^7re"^"^ { (r^ + l)x + 4r} 

+2t'^ {n - 2Si(a;)} {4cosj; + (r^ + l)2:sina::} 

-|-4T^Ci(a;) |4sin2: — (t^ -I- l)xcosa;} 



(B2) 



for n — 2, where a; = (1 -f z)krs. Note that mnfw(cvir) in 
the denominator originates from the division of the density 
profile by Mvir; if instead the profile is normalised by the 
total mass Aftot, ?Tinfw(cvir) should be replaced with mtot- 
The functions P{x) and Q{x) are defined as 



P{x) — sinh(x)Chi(x) — cosh(a;)Shi(x), 
dP 



Q{x) 



dx 



cosh(a;)Chi(a:) — sinh(a::)Shi(a::) 



(B3) 
(B4) 



For reference, the sine and cosine integrals and the hyper- 
bolic sine and cosine integrals are given by 



Si(a;) 

Ci(a;) 
Shi(a;) 
Chi(a;) 



sint 



dt. 



= 7 + In a; -f 



cos t — 1 

t 



dt. 



sinht 
t 



dt. 



7 + In a; + 



cosh t — 1 
t 



dt. 



(B5) 
(B6) 
(B7) 
(B8) 



with 7 = 0.57721566 . . . being the Euler-Mascheroni con- 
stant. While P{x) and Q{x) can easily be computed numer- 
ically (except at large x where the round-off error domi- 
nates), we provide fitting formulae for them: 



Pflt(a:) 



+ 



6a: 



c+(a:-d)2 
-|-a;(7 -I- In a; — 1) ( — 



X 



a;4 + a4 
4 \f 



(B9) 



with (a, 6, c, d, e, /)=(1.5652, 3.38723, 6.34891, 0.817677, 
-0.0895584, 0.877375), and 

9 



Qfit(a;) = 



+ 



bx" 



c+{x~ d)* 



i-y + lnx) ( 1 + ^J - 



x-^ + a-* 



(BIO) 



with (a, b, c, d, e, /, 3) = (2.26901, -2839.04, 265.511, 
-1.12459, -2.90136, 1.86475, 1.52197). These fitting formu- 
lae are obtained simply by connecting series expansions at 
a: = and 00. We find that Pfit(x) is accurate at better than 
0.5% for any positive values of x. The function Qmix) is also 
accurate at better than 0.5%, except at around x ~ 0.88 
where Q{x) has a root (i.e., Q{x) « 0). 
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